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Abstract 

The intrinsic volumes induced by a stationary Poisson fc-fiat process inside a com- 
pact and convex sampling window are considered. Using techniques from stochastic 
analysis, more precisely calculus with multiple stochastic integrals and a Wiener- 
Ito chaos expansion of Poisson functionals, all moments and cumulants, exact and 
asymptotic, are calculated in terms of a family of integral-geometric functionals. 
Moreover, univariate central limit theorems as well as Berry-Esseen-type inequalities 
are shown and a multidimensional limit theorem is concluded. 
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1 Introduction 

Random subdivisions of the d- dimensional Euclidean space (with d > 2) induced by a 
collection of random hyperplanes are one of the most classical and well studied models for 
random tessellations considered in stochastic geometry and spatial statistics. A class that 
has attracted particular interest in the literature is that of Poisson hyperplanes, which 
was introduced by Miles [9] and Matheron [8] and was studied by many authors (see 
[15] and [19] and the references therein). Also the more general class of k- flats (0 < 
k < d — 1), in particular Poisson Axfiats, was considered in the literature [1, 4, 15, 17, 18]. 
However, only a very few explicit and non-asymptotic results are available for these models. 
Main contributions in this direction are due to Heinrich et al. [3] and Heinrich [2] in the 
hyperplane case. The aim of the present paper is to exploit a new technique resting upon 
a Wiener-Ito-chaos decomposition and Fock space representation of Poisson functionals 
[7] and to show how it can successfully be applied to get new and deep insight into the 
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structure of Poisson fc-flats in M. d . We are able to extend the known results from the 
literature in two directions. On the one hand and in contrast to [2, 3] we will not only 
deal with hyperplanes, but with more general /c-flat processes and consider certain lower- 
order intrinsic volumes, whereas [2, 3] only deal with /c-volume and Euler-characteristic, 
two special intrinsic volumes, which are included in our discussion. On the other hand, 
we can derive elegant exact and non- asymptotic formulae for all moments and cumulants 
- completely characterizing thereby also distributions - in terms of a class of integral- 
geometric functional. The precise statements of our results are the content of the next 
section. 

As already signalled above, our results rest upon a Wiener-Ito chaos expansion of Poisson 
functionals. This is a very powerful and advanced technique from stochastic analysis that 
has not been used so far for problems in stochastic geometry. To keep the paper self- 
contained and for easier reference we review in Section 3 some of the most important 
results related to that theory, which are needed for our proofs. They are the content of 
Section 4 and the key to them is a product formula for multiple stochastic integrals. Such 
a formula is available in the literature in great generality, see [13, 20]. However, for our 
purposes a special and easily understandable particular case is sufficient. For this reason 
and in order to keep the argument simple, elegant and transparent, we circumvent the use 
of such general results and derive a recursion formula for the product of a special class of 
multiple Wiener-Ito integrals with respect to Poisson processes, which is much easier to 
handle than the classical product or diagram formula. 

We would like to point out that such an elegant and non-technical way is no more possible 
when dealing with lower- dimensional intersection processes and their induced intrinsic 
volumes. Their study is for this reason postponed to a future separate paper. 

2 Statement of results 

In this paper, we study stationary Poisson A;- flat processes r/ k in R d (0 < k < al — 1) with 
intensity measure r k K k) where r k e (0, oo), that is to say r\ k is a stochastically translation 
invariant point process on the space S k of fc-flats (this is the space of /c-dimensional affine 
planes) in M. d and the mean A;- volume EVolfc(% DT) (here and in the sequel Vol^ stands for 
the /c-dimensional volume measure) generated in a test set T C M. d of unit volume equals r k 
(as usual, we identify a point process with the random closed set induced by it). Here and 
in the sequel, A k stands for a unit density translation invariant measure on S k , by which we 
mean a measure of the product form £®TZ k , where £ stands for the Lebesgue measure on R 
and TZ k for a probability measure on the space C k of /c-dimensional linear subspaces of M. d . 
If lZ k is the uniform distribution on C k with probability element dL k , then rj k corresponds to 
a stationary and isotropic (this is stochastically rotation invariant) Poisson /c-flat process. 
The focus of our considerations is on the intrinsic volumes Vj^ k (W) = Vj(r] k nW), < j < k, 
i] k induces in a rf-dimensional compact convex body W C ~R d with interior points. For a 
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convex set K C M. d they are defined by 

V d (K) = Vo\ d (K), Vj(K) = / Vo\ 3 (K\L 3 )dL 3 , 0<j<d-l 

K 3 K d~j JCj 

where K\Lj stands for the orthogonal projection of K onto Lj and where Ki stands for the 
volume of the i- dimensional unit ball. If the affine hull of K C W 1 is fc-dimensional, Vk(K) 
is the /c-volume, 2V k -i(K) the (A; — 1) -dimensional surface area, ^^ Vi(K) the mean width 
and Vo(K) = 1 the Euler-characteristic of K. By additive extension, the intrinsic volumes 
are also well defined for locally finite unions of convex sets, cf. [15]. In this sense V kjk (W) 
is just the /c-volume or Vo tk (W) the number of fc-flats of r] k in W. We will also consider a 
sequence W e = gW of dilated windows in which rj k is observed. 

Under our invariance assumption, it is readily seen that EV kik (W) = r k Vo\d(W). In sharp 
contrast to the mean values, higher-order moments will depend on the precise shape of the 
window W via integral-geometric functionals A(W,j, k,m) defined by 

A(WJ,k,m):= [ Vp(WnE k )A k (dE k ), m = 1, 2, 3, . . . , < j < k, 

J[W\k 

where [W] k C £ k denotes the set of /c-flats hitting W. Let us denote by ii m (Vj ;k (W)) the 
m-th moment of V jik (W) - EV jtk (W), i.e. fi m (V jtk (W)) := E(V j:k (W) -EV jtk (W)) m and by 
lm(Vj,k(W)) the cumulant of order m (recall that the cumulants of a random variable are 
the coefficients of the series expansion of the logarithm of its moment generating function, 
see e.g. [16, p. 290] or (5) below), in particular YVj t k(W) = ^(Vj^W)). We also consider 
the normalized intrinsic volumes V* k (W) defined by 

To neatly formulate our results, denote by II({1, . . . , m}) the set of partitions of {1, . . . , m} 
in disjoint sets with at least two elements and by U m the multiset 

II m = {{mi, . . . , m r } : (A)I=i e H({1, . . . , m}),mi = \ Ai\ for i = 1, . . . , r} . 

The cardinality of {mi, ... , m r } in the multiset Ii m is given by the multinomial coefficient 
(my 

\m\...m T l 

Theorem 1 Consider a stationary Poisson k-flat process in M d , < k < d — 1 as above, 
(a) For a fixed compact convex body W C M d with interior points we have fJ>i(Vjk(W)) = 

/U^W) = E ^ II Mw,j,k, mi ) 

vren m m ; e7r 

and 

-y m (V j>k (W))=T k A(W,j,k,m) 

for m > 2 and explicit expressions for A(W,j, k,m) are provided in Lemma 1 below 
for the case of a Euclidean ball. 
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(b) For the growing sequence W e = gW the following holds: 



lim 



^ m {V hk {W e )) l(m-l)\\(r k A(W,j,k,2))^ : m even 



^jm+i^ 1 } ' \(^){m-A)\W k A{W,j,k,3){T k A{W,j,k,2)) m 2 1 : m odd 
for m > 2 with the double factorial kU = k{k — 2){k — 4) . . . and 

(c) For the normalized intrinsic volumes V* k (W Q ) we have 

i- /t /tt7" \ \ J(m-l)!! : m even 
hm Hm{V* k {W e )) = \\ 

: m odd 



for m > 2 and 

lim lm{V* k {W e )) 



Q— >00 




In general, the integral-geometric functionals A(W, j, k,m) are very difficult to evaluate 
explicitly. However, for a stationary and isotropic Poisson /c-flat process r\ k and the d- 
dimensional ball B d Q with radius g > we have the following closed formula: 

Lemma 1 For an isotropic Poisson k-fiat process, W = B d , 0<j<k<d — 1 and 
m — 1, 2, 3, . . . we /iawe 



A(B d Q ,j,k,m) = 



k\ n k 



J J K k . 



m 

K k m 



Thus, in the isotropic case and for W = B d := Bf, the moments and cumulants of Vj >k (W) — 
EVj 5 fc(W) can explicitly be evaluated. It is also worth pointing out that in the isotropic 
set-up the special functional A(W, j, k, 2) may be interpreted as the order j + 1 chord-power 
integral of W. Indeed, a repeated application of the affine Blaschke-Petkantschin formula 
from integral geometry shows the equality 

A(WJ,k,2) = ^- f Voli+^WnLjdL!, 
J + 1 Jd 

where C\ is the space of lines in M. d with invariant measure dL, see [15, Eq. (8.57)]. 
It is easily verified that the moments fJ- m {Vj tk (W)) satisfy the assumption of the well known 
Carleman criterion for the uniqueness of the moment problem [16, p. 296]. Thus, the entire 
distribution of the centred intrinsic volumes - and hence the intrinsic volumes itself - is 
determined by the sequence f^ m {Vj^ k (W)). Further, the asymptotics for the moments show 
how the leading term of their expansion is influenced by the shape of the window W. 
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The knowledge of the distribution of Vj t k(W) allows alternative proofs of a central limit 
theorem for V* k (W g ), a topic which was previously discussed in [2, 3, 11] in the special cases 
k — d— 1 and j = or j = k. It is readily seen that the well known method of moments can 
be applied and also the celebrated cumulant method. Alternatively, the general estimate 
for the Wasserstein distance recently obtained in [12] by means of Stein's method combined 
with Malliavin calculus for Poisson point processes can be used to conclude the central limit 
theorem. Moreover, we are in the position to prove a multivariate central limit theorem 
for the vector (V^* fc (W e ), . . . , V kk (W Q )) of intrinsic volumes of a stationary Poisson fc-flat 
process, which is new. 

Theorem 2 Let rjk be a stationary k-flat process in M d , < k < d — 1. 

(a) As g — >■ oo ; V* k {W e ) converges in distribution to a standard normal random variable 
Af and, moreover, 



sup < t) - nu < oi < ««|r;V 
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(b) The vector (V^ k (W g ), . . . , V k k {W Q )) converges in distribution to a (k+1)- dimensional 
Gaussian random vector with mean zero and covariance matrix C = (cjj)* - =0 with 
Cij given by 



C i,j 



•V Vi( - W n E k) y j(W n E k )A k (dE k ) 



V^W V ^ W n E *)MdEk) f [wh Vf{W n E k )A k (dE k ) ' 
In the particular case W = B d we have 



r(i + i±*) / r(j + i)r(l+j) 
r (1 + i?) V <yi 



Remark 1 Beside the multivariate normal convergence of part (b) in the previous theorem, 
also a rate of convergence for certain multivariate probability metrics can be concluded from 
[14]- However, their definitions are quite involved and the corresponding results are omitted 
for this reason. 



3 Wiener-Ito Chaos Expansions 

This section is devoted to a brief introduction to Wiener-Ito chaos expansions. Let 
(X,B(X), X) be a standard Borel space, which means that there is a bijection between 
X and a Borel subset of [0, 1], which is measurable in both directions, see [5]. Moreover, 
we assume that the measure A is a-finite and has no atoms. Let r\ be a Poisson point 
process with intensity measure A on an underlying probability space (Q, J- ', P) . By we 
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denote the distribution of r] on the space N(X) of all integer-valued cr-finite measures on 
X. The space N(X) is equipped with the smallest u-field such that all maps v h-> v(A) 
with A G B(X) and v G iV(A') are measurable. 

We present now a definition for multiple Wiener-Ito integrals, following thereby [20]. One 
starts with simple functions and extends the definition to arbitrary functions in L 2 (X n ), 
n G N, the space of square-integrable symmetric functions with norm || • || n and inner 
product (•, •)„. A function / G L 2 (X n ) is called simple, if 

1. / is symmetric, 

2. / is constant on a finite number of Cartesian products B ± x . . . x B n and vanishes 
elsewhere, 

3. / vanishes on diagonals, that means f(x±, . . . , x n ) = if x m = xi for some m ^ I. 

Let Ll(X n ) be the space of all simple functions. For / G Ll(X n ), the multiple Wiener-Ito 
integral I n (fo) of ft with respect to the compensated Poisson point process 77 — A is defined 

as 

Wo) = / fod( V - A) = fo lX - XBn (v ~ A)(S0 . . . ( V - X)(B n ), 

where we sum over all Cartesian products and f^ x --- xB " [ s the constant value on such a 
set. A straight forward computation shows that 

E/ n (/o) 2 = n!||/o||^. (1) 

Thus, there is an isometry between Lq(X u ) and a subset of L 2 (P^). Furthermore, Lq(X u ) 
is dense in L 2 (X n ), whence for every / G L 2 (X n ) there is a sequence (fi)i<=n of simple 
functions with ft — > / in L 2 s (X n ). Because of the isometry (1), it is possible to define I n (f) 
as the limit of (I n {ft))ie® in L 2 (F V ). 

Remark 2 The fact that Ll(X n ) is dense in L 2 (X n ) depends on the assumed topological 
structure of the space X and the fact that A is non-atomic. For a definition without these 
requirements we refer to [7]. 

It follows directly from the definition that multiple Wiener-Ito integrals have the properties 
summarized in the following 

Lemma 2 Let f G L 2 (X n ) and g G L 2 (X m ) with m,n>l and n^m. Then 

(a) EI n (f) = 0, (b) EI n (f) 2 = n\\\f\\l (c) EI n (f)I m (g) = 0. 

The second ingredient for Wiener-Ito chaos expansions is the difference operator D. Let 
F : N(X) 4 1 := RU {i 00 } be measurable and x G X. In this case, the difference 
operator (or add-one-cost operator) D x is defined by 

D x F(rj) :=F( V + 5 X )-F( V ), 
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where 5 X stands for the unit mass Dirac measure concentrated at x. It is well known that 
the difference operator obeys the product formula 

D Z (FG) = (D Z F)G + F(D Z G) + (D Z F)(D Z G), (2) 

see Lemma 6.1 in [10]. The iterated difference operators D™ Xn with xi, . . . ,x n G X can 
now be introduced by 

D° = 1, D 1 := D, I)';. ,., !■■{!,) := D^D^^r,). 

Let us further define the functions T n F : X n h-> R via 

T F = EF(r)), T n F( Xl , ...,x n ) = ED^ Xn F(i]), n > 1, 

which are called the kernels of F in the following. Because of the symmetry of the iterated 
difference operators, the kernels are symmetric, too. It was shown in [7, Thm. 1.3] that 
any F G L 2 (P^) can be decomposed uniquely as a sum of random variables that are the 
multiple Wiener-Ito integrals of these kernels: 

Proposition 1 For F G L 2 {F V ) we have T n F G L 2 s (X n ),n> I, and 

F(,v) = E ^In(T n F) (3) 

n=0 

and the series converges in L 2 (P). Moreover, the kernels T n F are the X-almost everywhere 
unique functions g n G L 2 (X n ), n G N, satisfying F(r)) = Y^=i ^tjn(gn)- Furthermore, 

Cov(F, G) = J2 —(T n F, T n G) n (4) 

n=l 

for F,G<E L 2 (P V ). 

The space of sequences (f n )neN with f n G L 2 (X n ) and Yl^=i ^\\fn\\n < c>o is often called 
the Fock space in the literature. Note that by Proposition 1, we have an isometry between 
L 2 (P V ) and the Fock space. The identity (3) is called Wiener-Ito chaos expansion or Fock 
space representation of F. 



4 Proofs 

Recall that by r]k we denote a stationary Poisson fc-flat process with /c-volume density 
and that for a ci-dimensional compact convex body W, Vj t t(W) = Vj(r)k fl W) stands for 
the j-th intrinsic volume induced by rjk in W. Moreover, fJ> m (Vj,k(W)) ist the m-th centred 
moment E(V jjk (W) -EV j<k (W)) m and l m (V jik (W)) is the order m cumulant of V j>k (W), i.e. 

l 0g Ee^^W = £ ^l m (V^(W)), t G R. (5) 

777-. 
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By V* k (W) we denote the centred and normalized version (V j:k (W)-EV jik (W))/ y/VV j: k {W). 
Recall further that by W e we mean the family of dilated and growing windows gW C M d , 

Q > 0. 

In our exact moment formula, a crucial role is played by products of functionals A(W, j, k, m). 
To derive an asymptotic expression it will be necessary to identify the leading £>-term. To 
do so, the following observation will be helpful. 

Proposition 2 The functional A(W,j, k, m) is homogeneous of degree jm + d — k, i.e. 

A(W e ,j,k,m) = g^ +d - k A(W,j,k,m), g > 0. 

Proof Noting that Afc([WJfc) = g d ~ k A k ([W] k ) for translation invariant measures A k on 
E k implies 

A(W e ,j,k,m) 

= [ V J m (WnE k )A k (dE k )= [ Vp{g{WnE k ))g d - k K k {dE k ) 
J[w e ] k J[W] k 

= gi m+d ' k ! VJ n (WnE k )A k (dE k ) = ( f m+d ' k A(W,j,k,m), 
J[w] k 

where we have used that Vj is homogeneous of degree j. □ 

Proof of Lemma 1 Recall that we take for A k the translation invariant measure on S k 
induced by the unique unit density Haar measure on C k . In view of Proposition 2 it is 
sufficient to consider the case g = 1, where we write B d instead of Bf. Observe that the 
intersection B d f]E k (p) is a /c-dimensional ball with radius a/1 — p 2 within the intersection 
A;-plane E k (p), where p G [—1,1] denotes the distance from E k {p) e £ k to the origin. A 
direct calculation shows now 

A(B d ,j,k,m) = j\ m (B d nE k (p))dp= j 1 V J m (B k /T - I )dp 
= Vp(B k ) j\\-p 2 ) k -?dp = 

To the last formula we can apply Legendre's duplication formula for the Gamma function, 
which directly proves our claim. □ 

The proof of Theorem 1 rests upon the following results related to Wiener-Ito chaos ex- 
pansions. 

Proposition 3 (a) The functional Vj ik (W) has the Wiener-Ito chaos expansion 

V J , k (W)=EV J , k (W) + h(f 3 , k ) (6) 
with f hk {E k ) = Vj(W n E k ), E k G £ k . 



Km \ 
2 ) 

'! + ¥) 
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(b) All moments of Ii{fj, k ) exist and satisfy 

m 1 / — 1 \ 

EI l (f), k ) m = T k ^( m . jA(W;j,fc,< + l)E/ 1 (/ J - fc ) m - 1 - < (7) 

/or m > 2 . 

fc^ TTie recursion (7) has the unique solution 

^Hh k ) m = Yl T k l II Mw,j,k, mi ), (8) 

/or m > 1 . 

Part (a) tells us that ^-^(W) has a rather simple Wiener-Ito chaos expansion because all 
kernels T n S vanish for n > 1. As a consequence, the central moments of Vj >k (W) are given 
by E7i(/ Ji fe) m . As already discussed in the introduction, such powers of multiple Wiener- 
Ito integrals can be computed by a so-called diagram formula, see [13, 20]. But instead of 
this technically complicated formula we use another approach and show that the moments 
must satisfy the recursion (7) and derive a solution for this equation. 

Proof of Proposition 3 

(a) By the definition of the difference operator and the additivity of the intrinsic volumes 
one has for E k e S k , 

D Ek v hk {w){rik) = v j (wnE)-J2v j (wnE) 
= VjiwnE,). 

Since D Ek Vj ;k (W)(r] k ) depends only on E k and is independent of the Poisson fc-flat process 
Vk, 

D n w (n) V 3 , k (W)( Vk ) = 

k " " ' k 

for E£\ . . . , E { k n) e S k and all n > 2. Thus 

T^W^Ek) = ED Ek V jtk (W) = Vj(W n E k ) =: f hk {E k ) 

and 

T n V J 4W)(E { k 1 \...,E^) = 

for e£\ . . . ,E ( k n) e S k and n > 2. By (3), the Wiener-Ito chaos expansion (6) follows 
immediately. 



9 



(b) It is obvious that all moment of Vj tk (W) are finite, because < Vj tk (W) < cwjT]k(W) 
with a constant cwj > depending on the geometry of the window W and on the parameter 
j. For the proof of the recursion we start with the claim that 



z=i ^ ' 



(9) 



for E k e S k and with fj tk (E k ) = Vj(W fl i^). This is shown by the product formula (2) for 
the difference operator and by induction. For m — 1, the statement is easy to verify For 
this reason and to prepare the induction step, we start with m = 2 and obtain 

D Ek (h(fj, k ) 2 ) = 2h{f hk )D Ek h{f hk ) + (D Ek h(fj,k)) 2 = 2/,,( /•:,!/,(/,,) + f hk {E k f 
and moreover 

D Ek {h{f hk ) m+1 ) 

= {D Ek h{kk))h{kk) m + Ihihk) + D Ek (hU)M D Ek (h(f,, k ) m ) 

m f \ 



m+1 

E 

j=l 

m+1 

E 

i=i 



m 
i- 1 



m + 1 



+ 



fj,k(E k ) 1 h(fj,k) 



i=i x 

m+1— i 



fj,k( E k) l h(f : 



j,k) 



\m+l — i 



It follows directly from the series expansion (9) that T n (Ii(fj^ k ) m ) — for n > m — 1. 
Indeed, each application of the difference operator decreases the power of the Wiener-Ito 
integral by one. Thus, after an m-fold application of D, we get on the right side of the 
equation. To proceed, observe that by Lemma 2 and the Wiener-Ito chaos expansion of 
h(fj,k) m combined with the expansion (9), we get 



m(kk) m 
= E[/ 1 (/ iiJfc ) m - 1 / 1 (/ iiJfc )]=E 



m— 1 



nh{nD Ek {w hk ) m ' l ]))w h k)] 

m—l 



n ill \ k '"" 5 fc ' 

})^E k )¥,[D Ek {hU),k) m - l )}rk^k{dE k ) 



[w] k 



E 

i=i 

m—l 

r k J2 



m — l 

i 



[W] k 



f(E k y +1 T k A k (dE k ))E[h(f j:k 



\m— i— ll 



i=l 



jA(W,j,fc,i + l)E/ 1 (/ i>fc ) m - i - 1 . 
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(c) The identity 

i=i ^ ' 7ren m _i_i mien 

= e r i nl n ^j.^) 

follows directly by a decomposition of the partitions of m objects in a set which includes 
the first object and the partitions of the remaining objects. This finally completes the 
proof of the proposition. □ 

Proof of Theorem 1 By (8), we have 

^ m {V hk {W)) = nV 3 AW)-^V hk (W)T = m(fi,k) m = E T P II A(W,j,k,mi), 

7ren m m ; e7r 

which proves the first part of (a). 

The exact expression for the cumulants follows by a straight forward induction from the 
expression for /i m (S(W)) and by using the general well known recursion formula 



m 1 / — l\ 

lm(Vj, k (W)) = VmiVjAW)) " E ( 7- 1 MViAW^rn-iiVjAW)). (10) 

i=l ^ ' 



By the relation (10) between /i m and 7 m it is obvious that ^liVj^iW)) = (J>i(Vj >k (W)) = 
and -f 2 (V jik (W)) = fi 2 (V j:k (W)) - = r k A(W, j, k, 2). Furthermore, for m > 3 it holds 

i=l ^ ' 

= E ^ n^^^ m ')-Eu7i) r ^^^ A; ^) >< --- 

7rGn m+ i m^e^ i=2 ^ ' 

•••x E rPHAiWJ^m). 

7ren m+ i_i mieiT 

The second sum includes all partitions of the first sum with two or more sets, such that 

l m+ i{V hk {W)) = r k A(W,j,k,m + l), 
completing the proof of part (a). 

The asymptotic expression for the moments can be derived as follows. By Proposition 2, 
A(W,j, k,mi) is homogeneous of degree jmi + d — k, thus p appears in the expansion for 

M 

p m (Vj tk (W e )) with power ^^[jmi + d — k] = jm + \n\(d — A;), where \n\ is the number 
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of elements of n G IT m . For even m this is maximal if \n\ = m/2, i.e. if each mi in the 
decomposition of m is 2. For odd m, the decomposition of m into a sum mi+. . ■+m| 7r | must 
contain a summand 3, say m\ = 3, to be maximal. The remaining summands % > 2, 
can then chosen to be equal to 2. The combinatorial factors (m — 1)!! and — 4)!!, 

respectively, simply take into account the number of partitions which realize this situation. 
The asymptotic formula for the cumulants can be seen from 

lm(V j>k (W e )) = r k A(W g ,j,k,m) = r k ^ m+d - k A(W,j,k,m), 

where we have used Theorem 1 (a) together with Proposition 2. The formulas for the 
normalized intrinsic volumes V* k {W e ) follow upon division by (yVj >k {W e ))^ , in which g 
appears with power jm + m ( d ~ k ^ . □ 

Proof of Theorem 2 As already discussed, there are several ways to see the normal 
convergence of V* k (W g ). The characteristic function of V* k (W p ) can easily be seen to 
converge to e~* 2 / 2 , the characteristic function of Af, by using part (c) of Theorem 1. On 
the other hand, part (c) of Theorem 1 also shows that the cumulants of V* k (W g ) tend 
to those of a standard normal random variable. Another proof can be based on a general 
estimate for the Wasserstein distance of V* k (W g ) and M from [12], which is easily applied in 
our setting since V* k {W) has a simple Wiener-Ito chaos expansion. The Berry-Esseen-type 
estimate is a direct consequence from the result in [6] , saying that the rate of convergence 
in Kolmogorov distance is up to a constant given by 

T * f ,*w nP u/j^ T k A(W e ,j,k,3) 



[ V*(W g n E k )A(dE k ) = T k A(W e ,j,k,3) ^ 
(Wi, k (W e ))* J[w e ] k (T k A(W e ,j,k,2))z 



d-k 



which is of order r k 2 gr^r by Proposition 2. 
We turn now to the multivariate central limit theorem, part (b) of Theorem 2. By the 
chaos expansions of V^ k {W Q ) and Vj >k (W e ) and the covariance formula (4) with n — 1, we 
obtain 

Cav(Vi, k (W e ), V hk {W e )) = (f itk , f j , k ) 1 = [ Vi{W p n E k )Vj(W p n E k )r k A k (dE k ) 

J[Wp]k 

and normalization and the homogeneity in g leads to the claimed expression for Cov(V* k (W g ), V* k (W e )). 
In Corllaries 3.4 and 4.3 in [14], multivariate central limit theorems for vectors of first-order 
Wiener-Ito integrals are presented as a consequence of more general results. The bounds 
for two different probability metrics, which both imply convergence in distribution, are up 
to a constant, depending on the covariance matrix of the multivariate Gaussian distribu- 
tion, given by the sum over integrals of type (11) for i — 0, . . . , d, which also leads to a rate 

3. d—k 

of convergence of order r k 2 g~^~ . This immediately implies the statement by noting that 
the covariances &ov(V* k (W e ),V* k (W e )) do not depend on the scaling parameter g. The 
case W = B d can handled similar as in the proof of Lemma 1. □ 
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